High-concurrence steady-state entanglement of two hole spins in a quantum dot 

molecular 
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Entanglement, a non-trivial phenomena manifested in composite quantum system, can be served 
as a new type of physical resource in the emerging technology of quantum information and quantum 
computation. However, a quantum entanglement is fragile to the environmental-induced decoher- 
ence. Here, we present a novel way to prepare a high-concurrence steady-state entanglement of 
two hole spins in a quantum dot molecular via optical pumping of trion levels. In this scheme, the 
spontaneous dispassion is used to induce and stabilize the entanglement with rapid rate. It is firstly 
shown that under certain conditions, two-qubit singlet state can be generated without requiring the 
state initialization. Then we study the effect of acoustic phonons and electron tunnelings on the 
scheme, and show that the concurrence of entangled state can be over 0.95 at temperature T = IK. 

PACS numbers: 



Introduction- Semiconductor technology toward quan- 
tum information science has opened up the possibility of 
constructing scalable quantum devices. As an attractive 
host for storing quantum information bit (qubit), elec- 
tron or hole spins in self-assembled quantum dot(QD), 
are most promising for their scalability, relatively ease of 
coherent manipulations [l| and strong robustness against 
relaxation 0, In the past few years, significant theo- 
retical and experimental works have been made towards 
controlling and entangling quantum dots. These experi- 
ments include the efficient state initialization Q and co- 
herent population trapping of single spinQ, the spin- 
readout [2| and single-spin Faraday/Kerr rotations for sin- 
gle quantum dot spinQ, Q, as well as the inter-dot cou- 
pling in double quantum dots molecules @, Theoreti- 
cally, several schemes for entangling quantum dots have 
also been proposed [^, . 

However, all these schemes are exclusively tailed for 
electron spins in self- assembled quantum dot. Due to 
the longer coherence time compared with electron spin, 
hole spin in quantum dot has been paid more and more 
attention. A hole-spin state, which is constructed from 
a p-type atomic wave function, has many favorable as- 
pects such as highly-suppressed hyperfine interaction 11 
and much smaller tunneling rate[12|, in comparison with 
electron spin. Recently both experiments have been re- 
ported for initializing single hole spin with high fidelity of 
0.99 lllf . and creating the coherent population trapping 
state[13[. 

In this letter, we present a scheme to generate high- 
fidelity steady state entanglement of two hole in a cou- 
pled QDM. Our scheme is based on spontaneous emission 
where the coupling between the QDM is dominated on 
Foster resonant. We show that our scheme is initial-state 
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FIG. 1: (Color Online). (a) Sketch of the two vertically stacked 
quantum dots in QDM with Voigt geometry magnetic field. 
Forster interaction Vf exists between such two resonant QDs. 
(b) and (c) Four level scheme illustrating the ground and ex- 
cited states of a single self-assembled QD in the Voigt Con- 
figuration. (d)The preparation process of entangled state. 



independent, and robust against decoherence and tunnel- 
ing effect. At T — 1 K, concurrence higher than 95% is 
possible in state-of-art technology. 

The model - The QDM system composes two identity 
vertical aligned QDs, where an external magnetic field in 
the Voigt geometry [J] is used to break the degenerate of 
the energy level[14j (FigfTl (a)). Due to the build-in strain 
that from lattice mismatch between QDs and the host 
materials, the degeneracy of the heavy hole and light hole 
at r point is broken, hence only the heavy hole is only 
taken into account. Initially, each QD is dopped with one 
hole, where the initial spin polarization is not important 
in our scheme. The energy level of QD is shown in Fig. 
[I] (b), the ground hole states are |1) = | ft) and |0) = 



2 



| 4)7 an d the excited trion states are \s) = | ffjl-t) and 
\ r ) = I 1HH-)' where ft (JJ-) and | (J,) denote a heavy hole 
and an electron with spins along (against) x direction. 
As the magnetic field is along x direction, there is an 
additional Zeeman splitting E B ^ = g^^psBx between 
ground (excited) states. Because of this splitting, four 
transitions between ground states and trion states can be 
independently addressed by polarization and frequency 
selection. Here we can choose H = <r + — cr_ and V = 
0+ + cr_ . 

Optical pumping protocol - In generally, the tunneling 
rate decreases exponentially when increasing the distance 
between two dots, while the Coulomb interaction, includ- 
ing both static dipole coupling V xx and Forster interac- 
tion VpHEE^! decreases generally cx d~ 3 . So, it is pos- 
sible to choose a property distance, where the Coulomb 
interaction dominant the interactions between the QDM. 
The Static dipole coupling V xx is the bi-trion energy shift, 
while the Forster interaction between the QDM is a kind 
of inter-dot interaction. 

Due to the Voigt geometry magnetic field, the en- 
ergy detuning between two H type transitions are Ah = 
E B +E B , while the energy detuning between V type tran- 
sitions are Ay = E B — E B . In our scheme, we drive a H 
polarized laser, and the detuning Ah enables the transi- 
tions Hi or H 2 to be independently addressed, which has 
been achieved in experiments with high fidility[|. In the 
following we choose the transition Hi, thus the energy 
levels of QD employed in our scheme can be described by 
three states: |0), |1) and \s) (FigQ](c)). The frequency 
and Rabi frequency with H polarized laser are supposed 
to be oji and f2 respectively. Additionally, since the direct 
excitation of the transition |0) <-> |1) is forbidden, Q m can 
be realized by employing a Raman transition with large 
detuning to an auxiliary excited state. The Hamiltonian 
of the QDM reads 

H = [(m i (s\e luJ,t + n m \0) l (l\ + H.c.)+uj\ S ) l ( S \} 

i=l,2 

+ V F {\ls)(sl\ + \Os)(sO\+H.c.) + V xx \ss)(ss\. (1) 

Since Hamiltonian Eq.pQ is of the symmetry forma- 
tion, it is convenient to introduce symmetric state \Sij) = 
-^(\ij) + \ji)) and anti-symmetric state |Ay) = -^=(\ji) — 
= 0,1, s). With the aid of Forster interac- 
tion, an energy shift is generated between symmetric ex- 
cited states \Sos), | Sis) and anti-symmetric excited states 
\A 0s ),\Ai s ). 

When the H polarized laser is driven to pump Hi tran- 
sition with detuning A = — Vp, the transitions |Soi) *-» 
\Sq s ) and |11) <-> I Sis) are resonant in the rotating frame. 
In the case fi, Q m <C |Vf| -C V xx , the populations on bi- 
trion and anti-symmetric single-trion states are nearly 
equal to zero, and can be eliminated adiabatically. Us- 
ing the symmetric and anti-symmetric notation we intro- 



duced above, the scheme is reduced to a 6-state system 
(Figfjjd)). The effective Hamiltonian can be written as 

H eS = v / 2r!|ll)(Si s |+r!|S , oi)(So.s| + r! m |So s )<5i s | 
+ v / 2ar l |00)(5 i| + V2a„|Soi)(ll| + /i.c. (2) 

Then we take the lifetime of trion states into account. 
The photon emission occurs via an decay of the state |s) 
into |1) with Ti or into |0) with To. The total sponta- 
neous rate is assumed to be T = Tq +Ti. Then we derive 
a master equation within a Markovian process, 

p = -i[H eS , p] + J2(Lptf - ±{tfLp} + ), (3) 

i 

where L x = v / r^(|00)(S .s| + ^\S i)(Si s \), L 2 = 
- v /TV2|v4oi)(Si s |, L 3 = VrT(|ll>(Si s | + ^|S i)(S 0s |), 

U = v ^V2|^oi)(S 0s |. 

Steady state entanglement -The basic preparation cy- 
cle in our scheme works as follows. The lasers with 
Rabi frequency v^^m produce the transitions 1 00) «-> 
|Soi) «-> |11). Attribute to the Forster interaction, en- 
ergy splitting generates between the symmetric and anti- 
symmetric states. By tuning the laser frequency to res- 
onant with the symmetric single-trion states, the transi- 
tions |Soi) <-> \Sos) and |11) <-> |Si s ) are created. Thus 
the lasers couple the three ground states |00), |5oi) and 
1 1 1) to the single-trion states |So s ) and |Si s ), and leave 
the entangled state \Aoi) = -^(|10) — 1 1) ) decoupled 
from them. On the other hand, spontaneous radiation 
performs from \ s) to |0) and |1), which is dissipation from 
single-trion states to subspace Mi = {|00), |Soi), |H)} 
and M 2 = {\Aoi)}. If the excited states decay into M 2 , 
the process is terminated when one trion is dissipated; if 
they decay into M±, it will go through the cycle again. 
Therefore, after a period of time, the trions are dissipated 
and the system will go into a steady state |Aoi), which 
is the maximum entangled state we require. 

For a potential experimental system, we can choose 
a QDM including two resonant QDs. The wave func- 
tions for electrons and holes in each QD is supposed to 

have a Gaussian form 4>e/h ~ expf — ^ x ti V 1 — kjt], where 
' '11 M ± 

the characteristic lengths l\\(x) = y/h/(uj^±^Tn e /f,). We 

assume that l» e = 4.4 nm, Ihu = 4 nm, l± = 1 

nm. The Forster interaction [17[ read as |Vp| = 

li!ed 3 (ri~) 2 ^(f)' wnere £ i s the dielectric constant, d 
is the distance between the QDM, I 2 = 2/(1/% + 1/1%) 
and a ~ Ul -> /2m^E [l8|. In the case l e ,lh ^> l z , 
F(x) = f£ J* dt-j^g exp[-f], where v = 2 j^f t iy In our 
scheme, we set the distance between two dots in QDM 
is 9.5nm and \a\ = 1.6 nm, We can estimate Vp = —0.2 
meV and V xx = 3 meV. If the x-direction magnetic field is 
B x = 1 T, and g-factors are g* h = -0.29 and g* e = -0.46, 
the Zeeman splitting of trion and ground states arc about 
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FIG. 2: (Color Online), (a) Concurrence of the entangled 
state as trion-phonon effect is taken into account. The relative 
parameters: tt = 20 ficV, tt m = 9 ^eV, T = 1.2 fieV. (b) The 
characteristic time To as a function of Rabi frequencies Q and 



E% = -27.78 peV and £^ = -17.94 peV. Thus the de- 
tuning between Hi and Pi is |As| = 45.72 peV. Due 
to the requirement of frequency selection, the Rabi fre- 
quency of laser fl should satisfy f2 <C |As|, and £1 = 20 
peV is suitable. 

FigHU a ) shows the entanglement of the two hole spins 
in steady state measured by concurrence C as a function 
of time[19|]. The Red solid line represents the case only 
with spontaneous radiation, which shows a near unity 
concurrence at T = K after a period of time. Here we 
suppose that Ti = T2 = T/2. The input density matrix 
we consider is given by pi = |(|00)(00| + |Soi)(Sui| + 
|^4oi)(A)i| + |11)(H|)- After a period of time, the output 
state ideally generates as 1$/) = |Aoi). If the evolution 
time t — ► oo, the density matrix p — > |Aoi)(Aoi|. It is 
illustrated in Figf2fb) that the character time To, which 
expresses that the time for achieving the steady statefiol]. 
depends on spontaneous radiation rate T and the Rabi 
frequency of pumping field f2. If the rate of spontaneous 
emission is fixed, the character time can be shorten as 



Rabi frequency of laser field increases, and finally satu- 
rates approximal to a value To ~ 10/T, which is ten times 
that of the lifetime of trion state. The optimal charac- 
teristic time is appeared when we tune the coupling fi m 
to satisfy fi m = 0.45£l(shown in Fig^b)). 

For a real QDM, the rate of the spontaneous radia- 
tion from trion state to electronic state is about 1.2 peV, 
therefore To has a value about 5.5 ns. 

The effect of phonon interaction - In Fig[TJ we have 
taken exciton as the auxiliary state, which is vulnerable 
by the vibrational modes of the surrounding phonons. 
The interaction between acoustic phonons and excitons 
may be mediated by deformation potential coupling and 
piezoelectric coupling. Thus the phonon coupling matrix 
element 12 ll is 



where M. 



3q,i = 

CO 



e*- d J[Mf,p e (q) - M£,p h (q)], 



1,3 



MA 



I/O 



(4) 



2p\q\Vc, 

Following the Marko 



p e/h (q) = Jd 3 r|</yJV<i r 
vian approximation 22], the master equation of the 
density matrix in the interaction picture with respect to 
H, may be reduced into a Lindblad form 



P 



= J("MNi + l)D[Pi]p + ^[i^p], (5) 



where D[P]p = PpP^ — ^{P^Pp}+ is the decay operator 
of phonon effect, and iVj = [exp(w i /A:BT) — 1] _1 . J(u>i) 
denotes the phonon spectral density, and there are two 
kinds of J(w,) in our model which can be written as 



J±(w) = J dtt(l ± sinc( 



where Gd(w) = 
- Ph f. 



))[G d (fv) + G p (fv)], (6) 



D hPh f, G p (w) = 



Here the piezoelectric coupling is P q 



t sin 9M p V9 + 7 cos 28 — 2 cos 4cp sin 2 6*, in which M p de- 
notes piezoelectric constant (23l|. Combining the phonon 
effect into the Master equation Eq. (|3|), we will get nu- 
merical results of the concurrence shown in Fig[2{b). As 
temperature increases, the concurrence of entangled state 
decreases. The parameters for the phonons is taken from 
Ref. . 

The effect of tunneling effect - The tunneling ef- 
fect is assumed to be much smaller than the Forstcr 
intcraction|24j in above discussion. This assumption can 
be technically accomplished by increasing the distance 
between double dots but simultaneously the Forster in- 
teraction is suppressed. It is demonstrated that the tun- 
neling effect of hole is less than that of electron by one or 
two orders of magnitude [ll] which guarantees the stabil- 
ity of the hole-included ground state. Thus, we primar- 
ily consider the influence of electron tunneling on our 
scheme. Using standard WKB method, the tunneling 
rate can be estimated to t e ~ — J 8VJu> exp[— ^P-^l . with 
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FIG. 3: (Color Online). Concurrence of the entangled state 
for different electron tunneling rates: (a) as a function of time 
at T = IK. (b) as a function of temperatures. 



lu = A^/Wl/d. If the case is d = 9.5 nm, V' e = 680 meV, 
and m e = 0.067, the tunneling rate is t e = 1.9 meV, 
which can not be omitted compared to the Forster cou- 
pling \V F \= 0.2 meV. 

The exciton we discuss above in trion \s) is intra-dot 
exciton(an electron-hole pair lies in the same dot), such 
as \e\h\) or \e\h\), where e\{h\) denotes that one elec- 
tron(hole) generates in the ith dot. Due to the electron 
tunneling effect, the trion might include inter-dot exci- 
ton (|eJ/4) or l e 2^i))' an d we denote this kind of trion 
as \t). Here, we do not consider the spontaneous radi- 
ation caused by inter-dot exciton. Therefore the effect 
of inter-exciton gives a new contribution to the initial 
Hamiltonian describing by Eq. ([1]) as 



H t = Y^{ u t\t)i{t\+t e (\s)i{t\+H.c)). (7) 

i=l,2 

We move the combined hamiltonian including tun- 
neling effect into an interaction picture with respect to 
(w + VF)(|s)j(s| + \t)i(t\), and proceed by transform- 

i=l,2 

ing the single exciton part of hamiltonian into a new 
basis as |V>i(3)> = cos 0\S o{1)x ) - sin0|£ O (i)t), 1-02(4)) = 



sin0|So(i)x) + cos^|5o(i)i}, where 9 



iarccot( 27-) 



and the detuning 6 = Vf + w — uif These two eigien 
states and [^3) are degenerated at energy E\ = 
|(— 5 + \J^t 2 e + 8 2 ), while |^ 2 ) an d \i>4) are degenerated 
at energy E% = i(— S — y 'it 2 + S 2 ). If we select the 
frequency of pumping laser as u>i = u> + Vf + E\, and 
guarantee the condition il <C \Ei — the effective 
Hamiltonian becomes as 



H cS = Sl m (V2\00)(S 01 \ + V2\S 01 ){ll\ + \S 0s )(S ls \) 
+ rLcos8(V2\ll)(iP 3 \ + \Soi){^i\) + h.c. (8) 



We consider the influence of both tunneling effect 
and phonon interaction and follow similar method as in 
Eq.(|l]). The concurrence of hole spin entangled state as 
a function of time at T = IK is shown in Fig[3] (a). We 
can find that the concurrence is beyond 95%, and the 
phonon-exciton process might be suppressed by decreas- 
ing electron tunneling rate. When the electron tunneling 
rate is slow as t e <C \u> — idt\, the Eq. ([5]) can be reduced 
to Eq. @. It indicates that the small electron tunnel- 
ing docs nothing more than an energy shift of the single 
exciton states whose effect can be offset by tuning the 
Rabi frequency of the external laser field. Fig[3] (b) illus- 
trates the concurrence of the steady state as a function 
of experimental temperature T. The high concurrence of 
stationary entangled state is much less influenced when 
decreasing the temperature, and the concurrence is more 
susceptive to temperature if electron tunneling rate in- 
creases. In our scheme, the electron tunneling has a value 
of t e ~ 2 meV, which is much smaller than the energy 
gap between inter-exciton and intra-exciton(~ 20 meV). 
Thus, the hole spin entangled state remains robust even 
if electron tunneling effect is considered. 

Conclution - To sum up, we have shown that a station- 
ary entangled state on spins with high concurrence can be 
prepared in a quantum dot molecular by technically de- 
signing the spontaneous dispassion processes. The hole 
spin for its small inter-dot tunneling rate is more suit- 
able to encode qubit compared with electron spin in our 
scheme. We also discuss the influence of phonon-exciton 
interaction and electron tunneling effect on the entan- 
gled state. For the real experiment with t e = 2 meV, 
the concurrence of entangled state is still over 95% at 
T = IK. 

Acknowledgments. This work was supported by 
National Fundamental Research Program (Grant No. 
2009CB929601), also by National Natural Science Foun- 
dation of China (Grant No. 10674128 and 60121503) 
and the Innovation Funds and "Hundreds of Talents" 
program of Chinese Academy of Sciences and Doctor 
Foundation of Education Ministry of China (Grant No. 
20060358043). 



5 



[1] Ronald Hanson, David D. Awschalom, Nature 453, 1043 
(2008) 

[2] J. M. Elzerman et al, Nature 430, 431 (2004) 

[3] D. Heiss et al, Phys. Rev. B 76, 241306(R) (2007) 

[4] C. Emary, Xiaodong Xu, et al, Phys. Rev. Lett 98, 

047401 (2007);Xiaodong Xu, Yanwen Wu, et al, Phys. 

Rev. Lett 99, 097401 (2007) 
[5] J. Berezovsky, M. H. Mikkelsen, et al, Science 320, 349 

(2008) 

[6] Mete Atatiire, Jan Dreiser, Antonio Badolato and Atac 

Imamoglu, Nature Physics 3, 101 (2007) 
[7] Lucio Robledo, Jeroen Elzerman, et al, Science320, 772 

(2008) 

[8] M. Bayer, P. Hawrylak, et al, Science291, 451 (2001) 
[9] S. K. Saikin, C. Emary, et al, Phys. Rev. B 78, 235314 
(2008) 

[10] A Kolli, B. W. Lovett, et al, Phys. Rev. Lett. 97, 250504 
(2006) 

[11] Brian D. Gerardot, et al, Nature 451, 06472 (2008) 
[12] C. Emary, L. J. Sham, Phys. Rev. B 75, 125317 (2007) 



[13] Daniel Brunner, Brian D. Gerardot, et al, Science 325, 
70 (2009) 

[14] H. A. Kramers, Proc. Acad. Amsterdam 33, 959 (1930) 
[15] DaeGwi Kim, Shinya Okahara, Masaaki Nakayama and 

Yong Gu Shim, Phys. Rev. B. 78, 153301(2008) 
[16] S. A. Crooker, et al., Phys. Rev. Lett. 89, 186802(2002) 
[17] Emil Rozbicki, Pawel Machnikowski, Phys. Rev. Lett 

100, 027401 (2008) 
[18] P.Y. Yu and M. Cardona, Fundamentals of Semiconduc- 

tors (Springer, Berlin, 2005), 3rd ed. [10] W. K. 
[19] Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[20] B. Michaelis, C. Emary, et al., Europhys. Lett. 73, 

677(2006) 

[21] B. Krummheuer, V. M. Axt, and T. Kuhn, Phys. Rev. 

B. 65, 195313(2002) 
[22] Erik M. Gauger, Ahsan Nazir, Simon C Benjamin, 

Thomas M Stace and Brendon W Lovett, New J. Phys. 

10, 073016 (2008) 
[23] T. Takagahara, Phys. Pev. B 60, 2638 (1999). 
[24] Ahsan Nazir, Brendon W. Lovett, Sean D. Barrett, John 

H. Reina, and G. Andrew D. Briggs, Phys. Rev. B. 71, 

045334(2005) 



